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Using the CCFR data for the Gross-Llewellyn Smith (GLS) sum rule, we extract the strong 
coupling constant via Borel resummation of the perturbative QCD calculation. The method in- 
corporates the correct nature of the first infrared renormalon singularity, and employs a conformal 
mapping to improve the convergence of the QCD perturbation expansion. The important twist-four 
contribution is calculated from resummation of the perturbation theory, which is based on the ansatz 
that the higher-twist contribution has a cut singularity only along the positive real axis on the com- 
plex coupling plane. Thus obtained, the strong coupling constant corresponding to the central GLS 
experimental value is in good agreement with the world average. 



I. INTRODUCTION 



OA ' Many low-energy Quantum Chromodynamic (QCD) observables, including the Gross-Llewellyn Smith (GLS) sum 
rule, with characteristic energy scale a few GeVs are analyzed in perturbative QCD within the framework of operator 
product expansion (OPE). In this scheme usually the most important is the perturbative contribution from the 
Wilson coefficient of the unit operator, and there are nonperturbative, power suppressed, higher-twist contributions. 
•- Generally, the higher order coefficients in the perturbative contribution grow rapidly due to the asymptotic nature 
of the perturbative expansion. The uncalculated higher order corrections are thus expected to be large, and this 
£S) ■ can cause a large uncertainty in data analysis that employs the unprocessed, finite order perturbative expansion. It 
is therefore important to properly handle the divergent perturbation expansion via resummation, since it can give 
a more accurate result with reduced theoretical uncertainty. Besides, resummation serves to give a well defined 
O |' meaning to the higher-twist contributions. Without a proper resummation of the perturbative part, the higher-twist 
I \ contributions are ambiguous 

An often used resummation technique is the Borel resummation. It has a sound theoretical basis since it is built 
on our understanding about the singularities in the Borel plane which cause the divergence of the perturbative 
• • . expansion. Its use generally improves the quality of data analysis, as can be seen from the reduced dependence on 
m ^ ' the renormalization scheme and scale, and from the reduced dependence on the uncertainty of the uncalculated next 
, higher order perturbative coefficient. 

At moderate values of the strong coupling a s (Q) at a few GeV, the Borel integral receives most of its value from 
the interval between the origin and the first infrared (IR) renormalon singularity, and just beyond it, in the Borel 
plane. Let us call this loosely defined interval, for convenience, the primary interval. In Borel resummation it is thus 
very important to describe the Borel transform, which determines the Borel integral, as accurately as possible in the 
primary interval using the calculated first terms of perturbation theory. 

For this purpose two steps can be taken: (1) an explicit incorporation of the first IR renormalon singularity in 
the Borel transform, and (2) use of an optimal conformal mapping. With the usual power expansion of the Borel 
transform about the origin the information on the renormalon singularity is lost. To remedy this, one may explicitly 
incorporate the first renormalon singularity by writing a Borel transform D(b), which behaves as 1/(1 — b/bo) " 
around the singularity at b = bo, as 

with R(b) = D(b)(l — b/bo) 1+IJ . The function i?(6)[] is by definition bounded and has softer singularity at the first 



1 This function was first introduced in || in Borel resummation and independently in in renormalon residue calculation. 



1 



IR renormalon. Hence we can expect that the Borel transform in the form of (Q), with R(b) perturbatively expanded 
about the origin, would give a better approximation than the direct expansion of D(b). 

The step (2) utilizes the information on the locations of the singularities in the Borel plane. Use of conformal 
mapping in Borel resummation has a long history 0, and its use in perturbative QCD was particularly emphasized in 
H . On the Borel plane there are IR renormalon singularities on the positive real axis and ultraviolet (UV) renormalons 
on the negative axis.|| By pushing the singularities away from the primary interval one can obtain a smoother R in 
the new primary interval on the mapped plane. This would render the perturbation of R in the mapped plane to 
converge better. Among the several mappings considered in the literature we find that the one proposed in JoJ is 
particularly suited when used in combination with the Borel transform in the form of Eq. (|l|) . This mapping moves 
the first IR renormalon singularity that defines the primary interval to a point within the unit circle and all other 
singularities to the unit circle. Since the effect of the first IR renormalon is softened by the step (1), we expect this 
mapping in our case to be better suited than, for example, a mapping that moves all singularities to the unit circle 




These techniques were applied to the hadronic tau decay width l| and to the hadronic contribution to the muon 
anomalous magnetic moment In this work we apply them to the GLS sum rule. The CCFR analysis Jl^,[l4| of 
the GLS sum rule was based on evaluation of the truncated perturbation series (TPS) in MS scheme. Aside from the 
inherent ambiguity of the higher-twist contributions, this method gives predictions which are not stable under the 
inclusion of an additional term (~ a A s ) in the TPS. As we shall see, these problems can be avoided with the use of 
Borel resummation. 

A crucial new element of our analysis comes with the calculation of the nonperturbative contribution. Aside from the 
perturbative part, an important contribution to the GLS sum rule comes from the nonperturbative, hadronic matrix 
element of the twist-four operator. Being nonperturbative, this contribution is usually fitted using the QCD sum 
rule calculation. Recently it was proposed by one of us fl5[l , motivated by an observation that the nonperturbative 
amplitudes in lower dimensional solvable models have a simple analyticity in complex coupling plane, that these 
higher-twist contributions can in principle be calculated from the Borel resummation of the perturbation series. The 
proposal was based on the conjecture that the higher-twist contributions have cuts only along the positive real axis in 
the complex coupling plane, which allows to relate the real part of the nonperturbative amplitude to its perturbatively 
calculable imaginary part. This scheme was shown to work well in model field theories. When applied to some of the 
solvable lower dimensional theories, it allowed the associated nonperturbative amplitudes to be accurately calculated 
from the first terms of the perturbation theory in the respective theories. 

From our analysis we obtain for the strong coupling parameter the central value a s (Mz) ~ 0.117. Compared to 
the corresponding CCFR central value a s (Mz) = 0.114, our value is closer to the world average a 8 [Mz) ~ 0.118. 
The main improvement comes from the correct incorporation of the renormalon singularity on the Borel amplitude 
and the calculation of the nonperturbative contribution. 

The paper is organized in the following way. In Sec. |lj we describe the resummation method, incorporating in it 
the known structure of the leading IR renormalon and the nonperturbative part, as well as the conformal mapping. 
Section OH contains the numerical application of the method to the GLS sum rule, leading to predictions for a s (Mz). 
In Sec. [f7[ we compare our predictions with those of other methods. In Sec. |v| we discuss some general features of 
the Borel resummation and the OPE approaches to understand precisely where the two methods deviate from each 
other. In Sec. we summarize our results and present conclusions. 

II. THE METHOD 

In this section we give an overview of our method used for the QCD analysis of the GLS sum rule. Its implementation 
in detail will be given in the following section. The QCD correction A(Q 2 ) to the GLS sum rule is defined by 

f 1 dxF» N (x,Q 2 ) = 3(l-A(Q 2 )), (2) 
Jo 

where Fg N is the non-singlet deep inelastic scattering (DIS) structure function in vN scattering. Here we shall ignore 
the target mass correction since it is irrelevant for our present discussion, but it will be included in the numerical 
analysis presented in the next section. 



2 There are also instanton-caused singularities, which can be safely ignored in our case. 
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We first begin by reviewing the old but important problem Jl],|| with the conventional QCD formulation of A(Q 2 ) 
in OPE framework, which is widely used in data analysis. This problem is not confined to the GLS sum rule, but 
generic to any perturbative OPE formulations. A(Q 2 ) in OPE up to twist-four operator is given by 



MQ 2 ) = W u (n y[ Q)) + W\inAQ))^^ , (>!) 



where a s (Q) is the strong coupling constant, and <0> is the reduced nucleonic matrix clement of the twist-four 
operator that was first derived in Ref. [|l6| 

= \^ af} G a a ^ , (P\O^P) spin avoragcd = 2 Ptl « O » . 

Here, A a are the usual Gell-Mann matrices and we used the notations of Ref. Jr7| |. Throughout the article we shall 
consider only the twist-four contribution as the nonperturbative effect, and ignore higher twist contributions since 
they are believed to be small. 

In conventional QCD analysis the Wilson coefficients W,-, (i = 0, 1), in (|J) are taken from the finite order, pertur- 
bative QCD calculation in an usual renormalization scheme, say MS scheme, and the reduced matrix element <C03> 
from data fitting or QCD sum rule calculations, etc. However, this scheme is, in principle, fundamentally flawed, 
since perturbatively the Wilson coefficients are not well-defined. In perturbative calculation of the Wilson coefficients 
the quantum fluctuations of all energy scale contribute, and in particular at large orders the contribution from the 
far infrared regimes, where perturbative QCD should fail, is large and gives rise to a same sign, factorially growing 
large order behavior [l8| . Thus, without some kind of resummation of the divergent perturbation series, the Wilson 
coefficients are not well-defined, and therefore neither is the OPE ||). As a consequence, no well-defined meaning 
that is independent of the definition of the Wilson coefficient Wo, can be assigned to the matrix element <cO^ |jl|-f3). 

One could in principle introduce an infrared cutoff /i 2 (<C Q 2 ) in the perturbative calculation of Wo, and regard 
the twist-four term to contain all the low momentum — k 2 = K 2 < fi 2 contributions |l9|j2C| ]. This would remove 
the infrared renormalon and make the perturbation expansion for Wo convergent, and the twist-four term to be 
fi 2 dependent pOfl , However, to do this in practice is impossible, because it is impossible to compute arbitrarily 
complicated Feynman amplitudes of arbitrarily high power, and in particular their small momentum contributions, 
as stressed by the authors of Ref. [^0) . 

The problem discussed so far is not of academic nature only, but has important practical implications. One might 
still think the OPE ||), with finite order perturbative Wilson coefficients, is a good approximation scheme, since at 
any rate the perturbative Wilson coefficients can be regarded as a good approximation at reasonably small values 
of the strong coupling constant. Actually, this would be the case, provided that the nonperturbative, higher-twist 
effect is far larger than the ignored next higher perturbative term in Wq. In practice, however, this condition is not 
supported by data analysis. For instance, if this were the case, we would expect little variation in the fitted values for 
the twist-four contribution over the order of perturbation in Wq. But the variation is not small at all. In the QCD sum 



rule calculation using the next-to-leading order (NLO) Wq 1 17|, the twist-four contribution was found not to be small, 



roughly equal to the perturbative correction at Q = 1 GeV. But, it was observed in Refs. 21 12^] that the twist-four 
contribution virtually vanishes when fitted with the next-to-next -to-leading order (NNLO) Wo against the parton 
distribution functions extracted from experiments, and most of the twist-four contribution extracted in the NLO fit 
can be accounted for by the perturbative NNLO contributions. This can be interpreted as a clear manifestation of the 
inherent ambiguity of the OPE approach. Moreover, this tendency of strong dependence of higher-twist contribution 
on the order of perturbation appears not to be special to the GLS sum rule, but generic. The recent new estimate 
of the gluon condensate J2^| , from fitting the vectorial spectral function of hadronic tau decay using NNLO Adler 
function, gives a small central value of only 1/3 that of the original QCD sum rule estimate |m| which uses the leading 
order perturbation. These examples strongly indicate that the inherent ambiguity of the OPE can have important 
consequences in practical applications, and demands a careful treatment. 

Borel resummation resolves this problem, which proceeds as follows [Q . The perturbation series for Wo 

oo 

W = J2^na(Q) n+1 [a(Q)=a s (Q)/Tr] , (4) 
o 

which is, being of same sign at large orders^, non-Borel resummable at physical, positive coupling a(Q). So it is 



3 Here we ignore, for the moment, the UV renormalons, which give rise to sign-alternating large order behaviors. Being Borel 
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first Borel resummed at negative a(Q) which yields a Borel resummed amplitude Ap(a(Q)). Then to obtain a Borel 
resummed physical amplitude one may analytically continue Ap to positive a(Q) in the complex coupling plane. This 
gives for a(Q) > |§ 



i /*oo±ie 

A P (a(Q) ±ie) = — / dbe- b ^° a ^ W (b) (5) 

PO J0±ie 



with 



W (b) = Y / ^(b/(3oT- (6) 

Inserted for normalization convenience, po is the one-loop coefficient of the QCD /3-function. Eq. (||) shows explicitly 
that in case of existing singularities on the line of integration, i.e., IR renormalons, the real part of the Borel integral 
Ap is the (generalization of the) Cauchy Principal Value. The Borel transform (^|) is believed to be convergent on the 
unit disk \b\ < 1, and is known to have a branch cut along the positive real axis beginning at b = 1. Near the branch 
cut it behaves as || 

W (b) = fT^0o + " (1 - + 0(1 - b)) + (Analytic part) (7) 

with 

" = (Jh/0o-T2)/0o, (8) 

where (3\ and 72 are respectively the two-loop coefficient of the QCD /3-function and the one-loop coefficient of the 
anomalous dimension of the twist-four operator appea ring in the OPE (^). In our case v is positive; for instance, 
v = 32/81 when three quark flavors (nf — 3) are active |16)|] Note that the Borel transform beyond the convergence 
disk of the series (|6) can be obtained by analytic continuation. The analytic part, which is analytic around the 
singularity, is not ca culable, but the residue C, a real number, is calculable perturbatively f3|,p5[. 

Because of the branch cut the Borel integrals in (|5|) develop imaginary parts beginning at b = 1. Since the QCD 
correction A(Q 2 ) must be real, clearly Ap alone cannot reproduce the true amplitude. There must be something else. 
Precisely at this point, the nonperturbative amplitude, denoted Anp(o(Q)), comes to the rescue, which cancels the 
imaginary parts of Ap , rendering the sum of the two to be real. This Anp can be shown to be directly related to the 
twist-four contribution in the OPE (||) |p6| , with its overall form governed by the associated RG equation, and may 
be regarded as a refinement of the latter, now free from the inherent ambiguity of the OPE. The QCD correction may 
now be written as 

A(Q 2 ) =A P (a± ie) + A NP (a ± ie) 

= Re[A P (a±ie)] +Re[A NP (a±ie)] . (9) 

Either sign can be taken because the result is independent of the sign chosen. 

What can we tell about the nonperturbative amplitude Anp? Its imaginary part at a positive coupling is certainly 
calculable from the perturbation theory because it is essentially the imaginary part of Ap, albeit with opposite sign, 
which is calculable in principle from the perturbation theory. However, as we see in (|), what we need is the real 
part, which is certainly not directly calculable. 

There is, however, an intriguing possibility of perturbative calculation of the real part. From the perspective of 
Borel resummation the sole reason for the introduction of the nonperturbative amplitude above was to cancel the 
imaginary parts arising from the analytic continuation of Ap to physical coupling in the complex coupling plane. 
With the Borel integral (||) and the singular Borel transform (Q), we can see that Ap(a(<5)), which by definition can 
have a singu larity only along the positive real axis in the complex a(Q) plane, has a branch cut of the form (see 



Appendix A) 



summable, the UV renormalons can in principle be treated separately, and do not affect our discussion in any essential way. 

4 According to Ref. Jl(|: 72 = (N c —l/N c )/3 = 8/9. Our convention for parameters /3j (and Cj =/3j//3o) is specified by Eq. ( pf] ) 
in the next Section. For n/=3, we have /3o = 9/4, /3i=4. Therefore, ^ = 32/81. 
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- C(-a(Q))- u e^V/SoatQ) [ l + G(a) j _ (10) 

The nonperturbative amplitude Anp should cancel the imaginary part arising from this at positive coupling. The 
simplest functional form for Anp that can achieve this purpose can be obtained by postulating that, as has Ap, the 
nonperturbative amplitude have a branch cut only along the positive real axis in the coupling plane. This then leads 
to the following conjecture for Anp fL5[ : 

A NP (a(Q)) = C(-a(Q))- v e-V/MQ) [ i + (a) } (11) 

This has a very important implication because it allows to relate the real part of the nonperturbative amplitude to 
the calculable imaginary part. In this paper we will adopt this conjecture, which appears very plausible, at least to 
us, and as will be mentioned shortly it is supported by some lower dimensional solvable models. Now with and 
( |ll[ ) we can write the QCD correction with only the calculable Ap as 

A P (Q 2 ) + A NP (Q 2 ) = {Re T cot(^Tr) Im} [A P (a(Q) ± ie)] . (12) 



This equation is the basis of our numerical analysis in the following section. In Appendix A we present, for reference, 
some explicit formulas leading to (|io|)-(|l2"|). 

The argument that led to the determination of the nonperturbative amplitude above did not rely on any special 
property of OPE, but only on the general property of Borel resummation of a same sign divergent perturbation theory. 
Thus it can be applied to any perturbation theory with a same sign large order behavior. Application of this scheme to 
the lower dimensional, solvable models such as the double well potential and the two-dimensional nonlinear a model 
in large N limit, allowed an accurate calculation of the associated nonperturbative amplitude using only the first 
terms of the perturbation expansion in the respective models fl5| ]. Moreover, interestingly, the numerical estimate 
of the gluon condensate from this scheme applied to the Adler function gives a value virtually identical to the new 
estimate p3[ . 

Given the QCD correction in the form ( |l2| ) with Ap (a ± ie) given by (||) , our aim is now to describe the Borel 
transform Wo as accurately as possible in the primary interval using the known first terms of the perturbation series. 
For this, following the outline in Introduction, we write the Borel integral (|j) as 



Ap( ttW )±i £ )4f ±,£ ie«-i (13) 



with R(b) now defined by 



R(b) = (1 - b) 1+v W (b) (14) 



and v given in (^|) {y = 32/81). This step is expected to greatly improve the perturbative description of the Borel 
transform in the primary interval because it implements the renormalon singularity correctly, and renders us to deal 
with a much softer singularity. The singularity of R(b) at b = 1 is a branch cut and thus softer than that of Wo(b). 

We can obtain further improvement by use of a conformal mapping that exploits the known locations of the 
singularities of Wq (b) . The latter is known to have renormalon singularities at non-zero integer values of b on the real 
axis ^j,^8|. To speed up the convergence of the perturbative expansion of R(b), we may push the singularities, save 
the unavoidable first IR renormalon, as far away from the origin as possible. This way we can reduce the influence 
of the renormalon singularities and make R(b(w)) smoother around the primary interval. One such a mapping we 
consider is M 



VTTb-^j~b/2 8w 

WlO) = , —;^=^= => b(w) = . (15 

V1 VT+b+^/T^b/2 V; (3w 2 -2w + 3) V ' 

This maps the first IR renormalon to w = 1/3 and all other singularities to the unit circle \w\ = 1. We expect this 
mapping combined with the implementation ( pr|) to provide an optimized environment for the Borel integral. In the 
mapped plane the Borel integral now assumes the form 

A P (a(Q) ± ie) = i/ ^e-^/fr-W) fg£g) (16) 
00 J c± dw (1 - 6(w)) i+iy 

where one of the integration contours C_ is shown in Fig. 1. Again, since the answer is independent of the sign 
chosen, either contour can be taken. In the next section we perform a numerical analysis of our implementation of 
the QCD correction: Eqs. (pi) and (|l|). 
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III. NUMERICAL ANALYSIS 



In this Section we will apply the method described in the previous Section to the Gross-Llewellyn Smith (GLS) sum 
rule, deducing values of the QCD coupling parameter aif s (Af§ ) from the GLS values extracted from experiments. In 
the first Subsection, we will present the resummed expression for the contributions of the three massless quark flavors. 
As a byproduct of the obtained expression, we will obtain an estimate of the next -to-next-to-next -to-leading (N 3 L) 
coefficient w 3 of the perturbative expansion. In the subsequent Subsection, we will include the effects of the massive 
fourth quark flavor (c-quark) to the GLS observable. In the last Subsection, the available measured GLS values will 
be confronted with our resummed expression and values of a%®{M%) will be extracted. 



A. GLS — the massless nj = 3 part 



The GLS quantity M 3 (Q 2 ) is the following integral (first moment) of the charged-current non-polarized DIS 
structure function F 3 = (F^ P + F^ n )/2 over the Bjorken parameter x: 



M 3 {Q 2 ) = \ [ dx F 3 (x;Q 2 )^ 
J Jo x 



1 + 2 1 



^m N x 



2\ I/ 2 ' 



(17) 



where £ = 2x/ (1 + \J\ + 4m 2 N x 2 /Q 2 ) is the Nachtmann variable p9| , and tun is the nucleon mass. The quantity M 3 
is the first Nachtmann moment of F 3 which absorbs all the kinematical power corrections (target mass corrections: 
TMC) ~ {m 2 N /Q 2 ) n . The quantity Q 2 = —q 2 is the virtuality of the exchanged gauge boson; Q 2 characterizes the 
typical process momenta of the observable M 3 (Q 2 ). When expanding the integrand in Eq. (|j]) in powers of m 2 N /Q 2 1 
we can rewrite 



M 3 {Q 2 



dxF 3 (x;Q 2 



_ 2 m 2 N x 2 
3 Q 2 



2 2 \ 2 

m N x 



Q 2 



o 



' N 



Q ( > 



The above GLS moment can be written in the form 

M 3 (Q 2 ) = 3(l-A(Q 2 )) , 
where the "canonical" quantity A(Q 2 ) has the following power expansion Wo [cf. Eq. (§)]: 

A(<3 2 ) i-> W (Q 2 ) = a(l + w x a + w 2 a 2 + w 3 a 3 H ) . 



(18) 



(19) 



(20) 



Here, a = a s (/i 2 ; C2, c 3 , . . .)/tt is the QCD coupling parameter with a given choice of the renormalization scale (RScl) 
jji 2 and the renormalization scheme (RSch) parameters Cj (j > 2). The evolution of a with the RScl is governed by 



the renormalization group equation (RGE) 

da(n 2 ;c 2 ,c 3 , . . .) 
d In ji 2 



-(3 a 2 (l + cia + c 2 a 2 + c 3 a 3 + • • ■) , 



(21) 



where the parameters Cj = Pj/Po (j ^ 2) characterize the choice of the RSch, and (3q and c± = /3i//?o are universal 
constants.M The evolution of parameter a with Cj's (j > 2) is governed by analogous differential equations, which 
follow from the Stevenson equation (see Appendix A of Ref. [p0[ ) . The next-to-leading (NL) coefficient w\ has been 
calculated in Ref. ||l|] , and the NNL coefficient w 2 in Ref. [|32X~At the specific RScl ^i 2 = Q 2 and MS RSch, and when 
the number of the active quark flavors is n^ = 3, these coefficients have the values^ 



(o) _ 

■1 v ' 



w 1 {n 2 = Q 2 ) = 3.58333 



(0) 



w^^Q 2 ;^) = 20.2153 + w { 2 0) (I. I.) , 



u4 0) (U.) = -1.23954 . 



(22) 
(23) 
(24) 



Po = (11 — 2n//3)/4, ci = (102 — 38n//3)/(16/3o), where n/ is the number of active quark flavors. 



The superscript '(0)' 



„<") 



denotes the special RScl-choice /i =Q and MS RSch. 



G 



In the NNL coefficient, the small "light-by-light" part was separated off. The "light by light" contribution should 
not be included in resummations of A(Q 2 ; P), as will be argued at the end of this Subsection. 

The Borel integral (|l6|), which will be together with (|lj) the basis for our resummation, is independent of the choice 
of the RScl /i 2 and the RSch (c2, C3, . . .) used in the integrand.^ Thus, we can rewrite it as (see also Fig. 0) 



A P (Q 2 ±te) 



,±i( 



1 , db(w) 
dx- 



dw 



exp 



b(w) 



/3oa(M 



R(b(w);v 2 /Q 2 ) 



(1 - o(to)) 113 / 81 



(25) 



At this stage we can, as a byproduct of formula (0), obtain an estimate of the yet unknown N 3 L coefficient W3 
appearing in the perturbative expansion (^0|) [p|,|33|]. If working with the RScl /1 2 = Q 2 and in the MS scheme, the 
expansion of R(b(w)) in powers of w gives us: 



R(b(w)) = 1 + 0.526749 to + 0.709369 w 2 

+ (-43.2574 + 0.277464«4°V 3 + 0{w 4 ) 

(0) 



(26) 



where we excluded the "light-by-light" part (|24|) of the NNL coefficient w 2 (|2^) . Looking at the coefficients appearing 
it is reasonable to expect that the coefficient R3 at w 3 is |J?3| ~ 1. If we assume R3 — 1± 1, we obtain a rather 



stringent estimate w 3 — 159.5 ±3.6. If we adopt a more cautious assumption \R 3 \ ~ 10 1 , we obtain w 3 fts 160 ±30. 
If we apply Pade approximant (PA) [1/1] r{w) to the expansion (Eq) and re-expand it back in powers of to up to to 3 , 



we obtain an estimate 



,(0) _ 



159.3. On the other hand, if applying the [l/2]fl(to) and demanding to po i c = 1 (i.e., 



b = +2, IR2), the prediction is = 158.5; if demanding w vo \ c — — 1 (he., b= — 1, UVi), the prediction is w 3 v> = 157.0. 
Very similar estimates are obtained if we do not apply the conformal transformation 6(to) (|l5|). Therefore, we will 



(0) 



adopt the following estimate for w 



(0) 



,(0) 



[= to 3 (/i 2 = Q 2 ;MS)] = 158 ±30 . 



(27) 



We emphasize that this estimate excludes the "light-by-light" contributions, which are assumed to be suppressed 
at the N 3 L order. The exclusion of the "light-by-light" contributions reduces the (NNL) perturbative expansion of 
the (non-polarized) GLS sum rule to that of the Bjorkcn polarized sum rule (BjPSR) f32]. It is interesting that the 
effective charge method (ECH) p3-p6| and the (TPS) principle of minimal sensitivity (PMS) pTl] predict w 



(0) 



based on the assumption that cf CH 



„PMS 



Further, an RScl- 



is somewhat related to the PA and PMS approaches, also predicts 



130 

and RSch-invariant method |38j that 
u 3 ~ 130. This is at the lower end of our 
new estimate (|27|). We thus conclude that the explicit (and exact) structure of the leading IR renormalon of the 
GLS-observable (BjPSR-observable) in the Borel plane, as given in Eq. (0), is responsible for the somewhat higher 

estimate of w 3 °^ in comparison with the ECH, PMS, and PA-related methods of resummation. 

We add here a few remarks on the question of the "light-by -light " contributions. The power expansion of the 



massless part of the GLS sum rule to N 3 L order 
is given by 



(see Rcfs. 



in the MS scheme and at the RScl /i" = Q- 



{— A C F )W {Q 2 



= a(-^C F ) + a 2 C F (^C F 



3 

128 



1241 
"576" 



23 
16 
11 
12 



7 n f) 



5437 



864 



C 3 )C 2 C A + (-^- + ^( 5 )C F C 2 A 



+(- 



133 



,3535 



1152 24^"/ + ('1728 



115 



+ oCa " ^)C F C A n f - —C F 



864 



d abc d abc n 1 

+n/ ^V— ( ~192 + 8 C3 ^ 



0(a 4 



(28) 



where the Casimir coefficients for QCD (N c 
the last term in (E8) is 



= 3) are Cf=4/3, Ca = 3, and the group-theoretical factor appearing in 



7 The location of the renormalon pole b = 1 and the power v = 32/81 (ph are independent of the choice of the RScl and RSch. 
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d abc d abc = _ ^ 

3 

This group-theoretical factor is not present in the calculation of the GLS sum rule up to the two-loop (~a 2 ) order, 
and it appears for the first time at the three-loop (^a 3 ) order. This term is called "light-by-light," it corresponds 
to diagrams with a new topology involving exchange of three gluons. In our calculation we will add this "light-by- 
light" contribution A;j.(<5 2 ), given by the last term in Eq. (p8|), as a separate term not included in our resummation 
approach. The reason for this is the following: Resummation approaches cannot be expected to predict (and resum) 
those higher order terms which are characterized by new higher order group-theoretical factors, when we have only 
one such term (~ a 3 ) explicitly available. We assume that such resummed "light-by-light" contributions are small, 
comparable to the quite small ~ a 3 "light-by-light" contribution in (^8[). Similar considerations can be found in 
Refs. [p9|j37[, in cases of various observables and beta functions. 



B. Inclusion of the massive quark flavor (c) contribution and nuclear corrections 

In the approximation of massless quarks, the calculation of the N 3 L (~ a 3 ) QCD correction to the GLS sum rule 
has been carried out in Ref. |32j. The inclusion of the heavy (c-) quark contributions is important at the precision 
level at which we are working. In addition, it is important in order to estimate the scales Q 2 where the (c-)quarks can 
be treated as massless or massive quarks. This point is important because it indicates which number of light flavors 
n/ should be used in the (resummed) "massless" part of the the perturbation scries. The calculation of the heavy 
flavor contribution to the GLS sum rule up to the second order (~a 2 ) was performed in Ref. JiO) and discussed in 
pl| . According to this approach: (a) the quarks it, d, s are massless and result in the dominant nf = 3 massless QCD 
contribution to the GLS sum r ule A(Q 2 ); (b) for Q 2 » 2-4 GeV 2 , the massive flavor is the c-quark and it contributes 
a relatively small correction to the aforementioned n./ = 3 massless contribution. The other heavy quark flavors are 
ignored due to the strong suppression by the mixing angles in the Cabibbo-Kobayashi-Maskawa matrix and by the 
small values of Q 2 . 

The heavy (c-)flavor correction contributions to A(Q 2 ) are 



A C (Q 2 ) 
A C (Q 2 ) (1) 



A« + , 

-^—~a(n 2 )^ 
3(1 + 4 



, M1+|) 



A C (Q 2 )( 2 > = -a(n 2 ) 



> 2 Cf Tf 



16 




(30) 
(31) 



105^ 45 M 



6740 1 137552 1 62528 1 



A 5 \ 105 



102656 1 \ , ,A + 1. 
^05-^ ln( — 



63 £ 315 C 2 105 £ 3 , 
494 1516 1 23024 1 298432 1 
"63" + ~2\TJ + 63 I 2 + 315 I 3 " 



20 1 a 



A + l, 
A- V 



(32) 



if)— ( 32 ) will not be included 



where £, = Q 2 /m 2 , A = ^Jl+Am 2 /Q 2 , Cf = ^/S, Tp — 1/2. The heavy flavor corrections 
in the resummation procedure for A, but will be added separately. 

In addition to the discussed heavy flavor contributions, there are also nuclear corrections to the GLS sum rule, due 
to the nuclear effects in the F3 structure functions. These effects were calculated in Ref. |4^] and were found for the 
iron target (used by the CCFR Collaboration) to be small: 



A Fo (Q 2 ) 



4 • 10~ 3 GeV 

Q 2 



(33) 



This contribution, at Q 2 — 2-3 GeV 2 , turns out to be in its magnitude smaller than the heavy flavor effects by a 
factor of 3-4, but larger than the "light-by -light" contribution by a factor of about 2. 
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C. Extraction of a s values 



The values of the GLS sum rule M 3 (Q 2 ), for various specific Q 2 , have been obtained from the experimentally 
extracted values of the structure function F 3 (x;Q 2 ), by the Fermilab CCFR Collaboration They have large 

systematic experimental uncertainties, primarily because of the uncertainties in the normalization of xF% and in the 
integration in the regions x <C 1 and x > 0.5. The systematic uncertainties are somewhat smaller, comparable to the 
statistical experimental uncertainties, only when the values of the exchanged boson virtuality are Q 2 = 2.00, 3.16 GeV 2 
(see their Table III). These values, including the target mass correction terms of (|T7|)-(|l8|), as well as the nuclear 
correction term (|3^), are[] 

Q 2 = 2.00GcV 2 : M 3 (Q 2 ) = 2.49 ± 0.13 <^ A(Q 2 ) = 0.168 ± 0.043 , (34) 
Q 2 = 3.16 GeV 2 : M 3 (Q 2 ) = 2.55 ± 0.12 & A{Q 2 ) = 0.149 ± 0.039 . (35) 

Here we added in quadrature the statistical and the systematic experimental uncertainties. In addition, the values 



A(Q ) in (|34|)-(|35|) have the aforementioned nuclear correction contribution (33) subtracted out; however, A(Q ) still 



includes the massless (perturbative and nonperturbative) contribution fll2|), the light-by-light contribution W2(ld.) a 3 
[cf. (fj)] and the c-quark contribution @-(||). 

For a given value of a s (Q 2 , MS), our expression for A(Q 2 ) = A P (Q 2 ) + A NP (Q 2 ) + A LL (Q 2 ) + A C (Q 2 ), given in the 
previous Subsections, still has the freedom (ambiguity) of the choice of the RScl /x 2 and the RSch. 

In the following analysis, we fix first the RSch to be MS, with the beta function being the PA f3(x) — [2/3] (x) based 
on the N 3 L TPS of the (3^(x). This PA choice of the beta function has reasonable behavior in the region of large 
x = a s /n and has been used in Refs. |) 11 3f| in the analyses of low energy QCD observables. We will comment later 
on how the results change when we use N^L TPS of the f3j-^(x), and when we change the RSch even more drastically. 

The RScl-dependence is yet another source of the theoretical ambiguity. In Figs. and|3| we show the dependence 
of the predicted values of A(Q 2 ) on the RScl j ji 2 a t the NNL and N 3 L level, at Q 2 "^ 2.00, 3.16 GeV 2 , respectively, 
for given fixed representative values of a s (Q 2 , MS). The NNL level means that we take for R(b(w)) in the Borel 



integrations in (25) and (h2j) the NNL TPS (quadratic polynomial in w), i.e., only the knowledge of w\ and W2 
coefficients is used. The N^L level means that R(b(w)) is the N 3 L TPS (cubic polynomial in w), i.e., we use, 
in addition, wf ' 1 — 158 according to the estimate (|27|). From Figs. |2|J| we see that the N 3 L expressions drastically 
improve the stability of the predictions under the change of the RScl. In the case Q 2 = 2.00 GeV 2 , the N 3 L values of the 
massless quantity Ap + np(Q 2 ) achieve minimal In /i 2 -sensitivity at /J. 2 /Q 2 ~ 3.3 : dA P+ ^ P (Q 2 ) /din [i 2 » —1.45 • 10~ 4 . 
At Q 2 = 3.16 GeV 2 , there is unfortunately no point of minimal RScl-sensitivity; the slope is negative and getting 
weaker when [i 2 increases. However, the RScl-sensitivity is very weak when fJ, 2 /Q 2 > 2.5, and almost stabilizes 
when [i 2 /Q 2 « 3.3 : dAp + Np(Q 2 )/dln/x 2 sa —1.55 • 10~ 3 . Therefore, we will take, for definiteness, the RScl choice 
fi 2 = 3.3 Q 2 in the case of both Q 2 values. Later on, we will comment on the ambiguity of our results when the RScl 
is varied. 

Having fixed the RSch (MS, with [2/3]^ (af) PA) and the RScl (^ 2 = 3.3 Q 2 ), our expressions for A(Q 2 ) = A(Q 2 ,P+ 
NP + Z.Z. + c) become unambiguous functions of the input value of a s (Q 2 ;MS). Adjusting simply the latter value so 
that the experimental values (Em])-([35|) are achieved, we obtain: 



Q 2 = 2.00 GeV 2 : a s (Q 2 ;MS) = 0.348±o.o9 3 (exp.) (36) 
Q 2 = 3.16 GeV 2 : a s (Q 2 ;MS) = 0.305^;^(exp.) (37) 

The central, upper and lower values here correspond to the pertaining experimental values of A(Q 2 ) in (|3^)-(|35|). In 
the case Q 2 = 2.00 GeV 2 , the upper bound for the coupling parameter cannot be obtained from the experimental 
data because the applied resummation method predicts values of A(Q 2 ) which are always lower than the presently 
allowed experimental upper bound: A(<3 2 ) mct h. < 0.196 [< A(Q 2 ) max oxp . = 0.168 + 0.043]. The situation in the case 
Q 2 = 3.16 GeV is somewhat similar, the experimental upper bound being now slightly below the maximal value 
allowed by the method. This situation is presented graphically in Figs. || and || which show A(Q 2 ) as function of 
a s (Q 2 ) as predicted by the applied method, for Q 2 = 2.00,3.16 GeV 2 cases, respectively. The central, upper and 
lower experimental bounds are included as straight dotted lines. 

We now return to the question of the uncertainties of our predictions under the variation of the RScl and RSch. If 
we vary the RScl parameter £ 2 = /i 2 /Q 2 around 3.3 across the interval 1.5 < £ 2 < 5, the predictions for a s (Q 2 ;MS) 



8 In the entries of Table III of Ref. Jl3[, the target mass corrections are included, but the nuclear corrections neglected. 
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vary by at most ±0.002. On the other hand, variation of the RSch leads to larger ambiguities. For example, if we 
repeat the calculation in the 't Hooft RSch [c% = C3 = 0, we take [2/3]p(x) for definiteness)|] keeping £ 2 = 3.3, the 
predictions for a s (Q 2 ;MS) change by 0.009 (Q 2 = 2.00 GeV 2 ) and 0.005 (Q 2 = 3.16 GeV 2 ). We will regard these 
values as characteristic values for the RSch-uncertainties of our results. The replacement of the MS [2/3]^(a;) by the 
N 3 L TPS MS beta function changes our predictions for a s (Q 2 ; MS) by only about 0.001. 

Yet another source of the theoretical uncertainty in our predictions may be the truncation in the (N 3 L) TPS 

R(b(w)). We regard the uncertainty ±30 in the estimated value of Eq. (p7f), as the major source of the 

truncation uncertainty. This changes the prediction for a s (Q 2 ;MS) by only ±0.001. 

T he ot her sources of theoretical uncertainty come from the massive (c-)quark contributions presented in Sub- 



sec. |IIB| due to the uncertainties m c = 1.25±0.10 GeV and sin(6> Ca bibbo) = \V cd \ = 0.224±0.016 |43|. The resulting 
uncertainties of the predictions for a s (Q 2 ; MS) are ±0.002 and ±0.002, respectively, when Q 2 = 2.00 GeV 2 , and ±0.001 
and ±0.001 when Q 2 = 3.16 GeV 2 . 

The final predictions for a s are presented in Table [|. In the Table, we presented the central predictions for 





Q 2 = 2.00 GeV 2 


=> M% 


Q 2 = 3.16 GeV 2 


=> M 2 




0.348 


0.1166 


0.305 


0.1167 


Sa s (>0, exp.) 


±? 


±? 


±0.132 


±0.0128 


Sa s (<0, exp.) 


-0.093 


-0.0115 


-0.074 


-0.0118 


Sa s (RSch) 


±0.009 


±0.0008 


±0.005 


±0.0006 


Sa g (RScl) 


±0.002 


±0.0002 


±0.002 


±0.0003 


Sa s (w 3 ) 


±0.001 


±0.0001 


±0.001 


±0.0002 


Sa a (m c ) 


±0.002 


±0.0002 


±0.001 


±0.0002 


5a s (sin(9 c ) 


±0.002 


±0.0003 


±0.001 


±0.0002 


8a s (evol.) 




±0.0003 




±0.0003 



TABLE I. Predictions for ct¥^(Q^) and af s (M 2 z ), extracted by the comparison of the results of the applied resummation method 
with the measured GLS values (p4[) , (p5|) . 

a s (<5 2 ;MS) when Q 2 = 2.00 and 3.16 GeV 2 , and the uncertainties of the predictions due to various sources. Further, 
we RGE-evolved these predictions to the canonical scale Af§ and included the results in Table Q. This RGE evolution 
was carried out by using the [2/3] Pade approximant of the four-loop MS TPS beta function, using the values of the 
four-loop coefficient c^{rif) [Q and the corresponding three-loop matching conditions for the flavor thresholds.^] 
When adding in quadrature the various theoretical uncertainties, the predictions of Table | can be summarized as 

Q 2 = 2.00GeV 2 : a s (Q 2 ; MS) = 0.348±o. 93( ex P-) ± 0.010(th.) , (38) 

a s (A/|;MS) = 0.1166^ 0115 (exp.) ± 0.0010(th.) ; (39) 

and 

Q 2 = 3.16GeV 2 : a s (Q 2 ; MS) = 0.305^;^ 2 (exp.) ± 0.006(th.) , (40) 

=> a s (Ml; MS) = 0.1167t^i 2 ^(exp.) ± 0.0008(th.) . (41) 

We see that the predictions of the applied method suggest that the experimental data on the GLS should be refined 
significantly in order to increase the predictive power for the QCD coupling parameter. 

Another observation, evident from Figs. ||-|[ is that the nonperturbative massless contributions (NP) to A(Q 2 ) are 
very significant, and negative. They have their origin, as explained in the previous Section, in the first IR renormalon 
singularity at 6=1, or equivalently, they correspond approximately to the d=2 massless power corrections (~ 1/Q 2 ). 
For Q 2 = 2.00, 3.16 GeV 2 , they lead to about 11%, 8%, decrease of the value of A(Q 2 ), respectively. This is to be 



9 For comparison, cf s = 4.471 cf s = 20.99, for n/ = 3. 

10 For details on the corresponding evolution uncertainties, we refer to Ref. j]]. They include the variation when the [2/3] 
Pade form of the beta function is replaced by the TPS form. 
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contrasted with the heavy (c-)quark contribution^ which is positive and leads to only about 3.6%, 3.1% increase of 
A(Q 2 ), respectively. If the c-quark contribution were not included, the central predicted values in (|3^)-(^l|) would 
change to a s (Q 2 ) = 0.367, 0.316 [a s (M 2 z ) = 0.1183, 0.1180] for Q 2 = 2.00, 3.16 GeV 2 , respectively. 

The small negative "light-by-light" contribution was separated from our resummation and then added as the term 

A;.;.(Q 2 ) ps wf\l.l.)a i (^ 2 ]C2, ■ ■ .)■ The "light-by-light" part of A LL (Q 2 ), by the special topology of the Feynman 
diagrams representing it, is a quasiobservable in the sense that it is RScl- and RSch-invariant. Thus, the coefficient 
w^\l.l.) ( |24| ) is the leading coefficient of that quasiobservable and is therefore unchanged under the changes of the 
RScl and RSch. The "light-by-light" part decreases A(Q 2 ) by only about 0.4% and 0.3%, for Q 2 = 2.00,3.16 GeV 2 , 
respectively. 



IV. COMPARISON WITH OTHER APPROACHES 

One may ask how crucial is the introduction of the conformal transformation (|l^) for obtaining the numerical 
predictions (|3§|)-([4i"|). If we repeat the same analysis, but this time without the conformal transformation, and 
keeping £ 2 = 3.3, we obtain 

Q 2 = 2.00 GeV 2 : a s (Q 2 ;MS) = 0.346±o.o 92 (exp.), a s (M|;MS) = 0.1163^. OU3 (exp.); (42) 
Q 2 = 3.16 GeV 2 : a s (Q 2 ;MS) = 0.304l°;^(exp.)), a s (M 2 ;MS) = 0.1165±°;°^(exp.). (43) 

These results are very close to the results (B8h (41 ). Thus, we see that the introduction of the conformal transformation 
( |l5| ) , which had the task of reducing the influence of the UV and the nonleading IR renormalons, does not influence 
significantly the predictions. Therefore, we can conclude that these renormalon singularities (at b = —1, —2, . . . and 
b = 2, 3, . . .) are in GLS numerically much less important than the leading IR renormalon singularity (at 6=1), even 
when no conformal transformation is introduced. 

We can ask how our predictions compare with those of other, alternative, OPE based methods which, in contrast 
with the method applied here, do not take into account explicitly the structure of the leading IR renormalon singularity 
in the Borel plane. 

One such an alternative method is the (TPS) PMS optimization of the perturbative contribution, which fixes the 
RScl and RSch in the TPS in a judicious manner J30| ] . Resummations of the GLS sum rule based on this method 
were theoretically and numerically investigated in 1992 by the authors of Ref. They were confronting the TPS 
results with the measured values, paying particular attention to the RScl- and RSch-dependence of the predicted 
values of a s (M z ; MS). For the nonperturbative massless (twist-four) d = 2 contribution, they employed the positive 
value as obtained in Ref. Jl7| [Anp(Q 2 ) ~ 0.1 GeV 2 /Q 2 ]. Further, the authors of Ref. jl6| accounted for the quark 
mass threshold effects (rij ) by introducing a judiciously chosen weighted average of n/ — 3, 4, 5.Q Furthermore, they 
used the GLS measured values available at that time [A(Q 2 = 3GeV 2 ) = 0.167± 0.027] which differ significantly from 
the presently available values (|35|). Their central value prediction was a s (M z ) = 0.115, which is lower than our central 
value predictions ( j39[ ) , (f4l|) . 

The CCFR group [ p-3|J14| carried out a numerical analysis similar to that of the authors of Rcf. , but with the 
newer, lower, experimental data @-@ for A(Q 2 ), and using in the TPS-part MS RSch (and RScl [i 2 = Q 2 ). 
Their central value is a s (3 GeV 2 ) = 0.278 (l4|0 and aJM z ) = 0.114, thus slightly lower than that of Ref. Q, and 
significantly lower than our central value predictions ([39]) and (fid]). The principal reason for this difference shall be 

discussed in the following section. Further, if they included in their method the N 3 LO term in the TPS, with as 
estimated in Eq. (|27j), the predicted value of a s (M z ) would decrease by about 0.002. 

Furthermore, the CCFR group mentioned that their central value increases to a s (3 GeV 2 ) « 0.305 |TJ] and a s (M z ) = 
0.118 |l||14|][^] when they set the twist-four (d= 2) contribution approximately equal to zero. Such higher-twist values 




of the heavy (c-)quark as 



11 The latter is, to a large degree, a d — 2 massive power correction ocm 2 /Q 2 , see Subsec. 

12 This is different from our approach, where we separately added the contributions 
corrections to the massless n/ = 3 GLS sum rule, as recently suggested in Ref. |kl|. 

13 The values of « s (3 GeV 2 ) obtained from our central values of Eqs. (Jq^j and (i^) are 0.309 and 0.310, respectively. 

1 We note that the RGE evolution of a s from Q to M z gives in our approach different results: a s (3 GeV ) =0.278(0.305) 
gives a s (Mz) =0.1124(0.1162) when using the three-loop or four-loop TPS /3-function, 0.1123(0.1161) when using the (four- 
loop) [2/3] PA /3-function - we use the corresponding two-loop and three-loop matching conditions for the flavor thresholds 
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for the GLS sum rule are suggested by the calculation by the authors of Ref. J48| based on an IR renormalon model 
with dispersive approach of Ref. Q , and also by the calculation by the authors of Ref. Q based on the bag model. 
It looks reasonable that the result with the aforementioned IR renormalon model method gives prediction rather close 
to our prediction, since our method also accounts for the IR renormalon contribution, although in a different manner. 
However, the calculations in Refs. jl8|,0 apparently do not give us a clear handle on how to treat the perturbative 
part, i.e., whether to take it as a LO, NLO, or NNLO TPS, or in any other form. This is in contrast with our method, 
where the perturbative and nonperturbative parts are clearly connected with each other. We discuss this aspect in 
more detail in the next section. 

At this point, we would like to point out that the PMS method has a signal casting doubts on its applicability in 
the discussed GLS cases- namely, the PMS RScl is very low in this case: /4ms w 0.203Q 2 . For Q 2 = 2.00, 3.16 GeV 2 , 
this implies the scales /j.pms ~ 0.64, 0.80 GeV, respectively, which may be too low for the application of perturbative 
approaches such as PMS. The same problem appears when applying the effective charge (ECH) method fp4|-^6[ to 
these GLS cases. 

The present world average for the QCD coupling parameter is a^(M 2 ) = 0.1173 ± 0.0020 by Ref. JsTJ and 
0.1184 ± 0.0031 by Ref. [^2| . Comparing this with our predictions ( |39| ) and (fill), we see that the method applied 
in the present paper gives us the central values which agree well with the present world average. We wish to point 
out that this agreement suggests that the method applied in the present paper for the nonperturbative massless 
correction to A(Q 2 ) is at least consistent with the experiments. If this correction were zero, or had the opposite sign, 
the obtained central prediction for a s (M|) would be at the lower edge or even outside the interval of the present 
world average. These considerations do not necessarily imply, but indicate, that the applied method gives the correct 
nonperturbative contributions. For a more definite statement in this respect, the experimental uncertainties in the 
GLS sum rule would have to be reduced significantly. 



V. BOREL RESUMMATION VERSUS OPE APPROACH 



In the discussions so far, we considered the amplitudes of Borcl rcsummation or OPE approach only at fixed values 
of Q 2 . Here we mean, for convenience, by the OPE approach the usual perturbative expansion plus a power suppressed 
term representing the twist-four contribution. In this section we consider them over a continuous range of Q 2 . This 
slight change of view will reveal the characteristic features of the two approaches, and enable us to better understand 
the cause of the significant difference in the extracted strong coupling constants seen in the previous section. 

We first note the remarkable stability of the Borel resummed amplitudes over the order of perturbation involved 
in their calculation. In Fig. fl (a)-(b), we plot, over the interval 1 < Q 2 < 10, in GeV 2 , the real part of the Borel 
resummed amplitudes for Ap(Q 2 ) using NLO, NNLO, and N 3 LO perturbations and the corresponding amplitudes of 
the ordinary TPS's in the OPE approach. In all of the plots in Fig. | the N 3 LO QCD /3-function was used in the 
running of the strong coupling. We also take the RG scale at /i 2 — Q 2 , since the RG scale dependence is sufficiently 
small for our present discussion. The aforementioned stability of the Borel resummed amplitudes becomes clear when 
the two figures are compared. Note the variation in the Borel resummed amplitudes is very small, whereas the TPS 
amplitudes have significant order dependence. While this stability is not completely unexpected, because the leading 
renormalon singularity is effectively softened by the use of the function R(b) in the Borel integration, the degree of 
the stability is still remarkable. This suggests that the renormalon-induced asymptotic behavior of the perturbative 
coefficients sets in quite early in perturbation, and that the use of R(b) and conformal mapping in Borel resummation is 
very efficient in handling the renormalon singularity. We also note in passing that the stability of the Borel resummed 
amplitudes for Ap(Q 2 ) + Anp(Q 2 ), which are not shown in the figure, is comparable to that of Ap(Q 2 ). 

In the previous section we have seen there is a significant difference between the Borel resummation and the OPE 
approach in the prediction of the strong coupling constant. While there are obvious differences in the two approaches, 
it was not clear what aspect of the Borel resummation is primarily responsible for the difference. Is it because of the 
perturbative part Ap(Q 2 ) or because of our specific implementation of the nonperturbative part Anp(Q 2 ), or both? 

To answer this question we plot in Fig. ^ (c) the Borel resummed Ap(Q 2 ) + Anp(Q 2 ) against the OPE amplitudes 
(NLO TPS) + 0.1/Q 2 , whose power term representing the twist-four contribution is from the sum rule calculation 



|ft5| ; other details given in Ref. ||. The CCFR Collaboration apparently uses an approximate three-loop RGE evolution [a 
truncated expansion in inverse powers of ln(^, 2 /Aj^g-)] and different matching conditions for the flavor thresholds, possibly of 
Ref. @, giving them the values 0.114(0.118). 
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|0 and (NNLO TPS) + 0.02/g 2 . The small power term in the latter was chosen for the amplitude to match 
the NLO OPE amplitude at large values of Q 2 in the plots. In the figure we first notice that the NLO and NNLO 
OPE amplitudes with a large difference in twist-four contribution agree reasonably well over the whole range of 
Q 2 considered. This implies that the higher-twist term in the NLO amplitude can be largely accounted for by the 
NNLO perturbative term, which is in qualitative agreement with the observation in Refs. pl|jj^] . This also shows 
that the use of the NLO sum rule calculation of the higher-twist contribution with a TPS of different order, which is 
not an uncommon practice, can be dangerous. Higher-twist contributions calculated at a given order of the leading 
perturbative contribution should never be used with a TPS of different order. 

On the other hand, the Borel resummed amplitude is in a reasonably good agreement with the OPE amplitudes 
at large Q 2 (> 4GeV 2 ), but deviates significantly at small momenta. Obviously, this deviation at small momenta 
explains the difference in the prediction of the strong coupling. Before we answer the origin of this deviation, we note 
that the good agreement of the two approaches at large momenta is a nontrivial result. Even though the amplitudes 
from the two approaches should agree at very high momenta (or at small couplings), since they have the same low 
order perturbations, the degree of agreement seen here is unlikely to be a random consequence. For instance, if our 
Anp(Q 2 ) had the wrong sign or were zero, then we would see a significant difference at high momenta (see the dotted 
plot for the wrong sign case) . Thus this good agreement of our Borel resummed amplitude with the OPE amplitudes 
at high momenta may be regarded as a partial support for our prescription of the nonperturbative part. 

Now back to the question of what causes the deviation at low momenta. It is not difficult to see Ap(Q 2 ) must be 
responsible for the deviation. The reason is as follows. Since the nonperturbative part Anp(Q 2 ) essentially behaves 
like a power suppressed termj^Jand Ap(Q 2 )+Anp(Q 2 ) is in agreement with the OPE amplitudes at large momenta, 
so should it be at low momenta, too, were Ap(Q 2 ) to behave like an OPE amplitude. Thus the primary cause of 
the difference in the predicted strong coupling constant must be that the Borel resummed Ap(Q 2 ) at low momenta 
cannot be parametrized in the form of an OPE amplitude. We can see this explicitly by looking at the plots in Fig. 
| (d) where, as an example, the NNLO Re[A P (Q 2 )] is plotted against an OPE amplitude (NNLO TPS) + 0.16/Q 2 . 
The power suppressed term in the latter was fixed so that the two amplitudes match at high momenta in the plots. 
Clearly, they deviate significantly at low momenta, with the Borel resummed growing more slowly than the OPE 
amplitude as the momentum is decreased. 

This suggests the OPE amplitude tends to overestimate at low momenta. We can easily see that this tendency 
arises from the bad functional form of the polynomial Borel transform (TPS of (^|)) around the renormalon singularity 
at b = 1. In Fig. fj] (a) the N 3 LO Borel transforms in the two approaches are plotted in the variable b. For b > 1 the 
Borel transform in our approach defined through R(b) in (|l4|) is complex, and its real part is plotted. It is obvious 
that the OPE Borel transform is badly broken around and beyond the renormalon singularity. When the coupling 
is small this is not a serious problem because the dominant contribution to the Borel integral ([&]) comes from the 
region close to the origin. However, as the coupling becomes larger the relevant integration region extends to the 
renormalon singularity, and beyond, and as we see in the plots, the OPE Borel transform can grossly overestimate at 
large couplings. The amplitudes obtained from these Borel transforms are plotted in Fig. ^ (b). As expected, at small 
momenta (large coupling) the OPE amplitude is larger than the Borel resummed. Note, on the other hand, at high 
momenta it is smaller than the latter. This is because the OPE Borel transform in the region 0< b< 1, from which 
the dominant contribution comes at small couplings, is smaller than the other one, which is a characteristic feature 
rendered automatically by the correct implementation of the renormalon singularity in the latter. This difference 
between the Borel resummed and the OPE amplitude at small couplings may be regarded as the resummation of the 
unaccounted higher order terms in the same sign asymptotic series. 

That the Borel resummed Ap(Q 2 ) cannot be approximated by an OPE amplitude of a TPS plus a power suppressed 
term representing the renormalon effect may appear contradictory to the common opinion which states otherwise. 
The latter opinion, which is based on the factorially growing large order behavior and the running coupling, would be 
true in a sense, provided the strong coupling were sufficiently small, and a TPS of large order (~ l/a s ) was used. In 
reality, however, the strong couplings at the low momenta we consider are not so small, and there is no guarantee that 
the Borel resummed Ap(Q 2 ) at those momenta can be parameterized as an OPE amplitude. The example here clearly 
shows that cannot be generally true. A Borel resummed amplitude can have a much more complex functional behavior 
than the sum of a TPS and a power term intended to account for the renormalon effect. This consideration suggests 
that several existing analyses of low energy QCD observables based on the OPE approach should be reexamined, since 



1 We note, however, there are some variations in the estimate of the twist-four contribution. The sum rule calculation of 
Ref. H predicts the (d = 2) power term 0.16 GeV 2 /Q 2 . 
16 A N p(Q 2 ) ~ a s {Q 2 y i2/l3 ° /Q 2 = a 3 (Q 2 ) 32/sl /Q 2 , in accordance with Eqs. (0)-© and the OPE calculation of Ref. (l(J. 
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the issues raised here are likely to be relevant there, too. 

To sum up this section, we have made two observations concerning the Borel resummation and the OPE approach. 
First, our method of calculation of the perturbative plus nonperturbative contribution in Borel resummation is 
consistent at larger Q 2 > 4 GeV 2 with the OPE approach using QCD sum rule calculation, and secondly, at low 
energies the OPE amplitude tends to overestimate, and the Borel resummed amplitude with a proper incorporation 
of the leading renormalon cannot be approximated by an OPE amplitude of a TPS with a power suppressed term. 
It is the second observation that directly accounts for the differences in the extracted strong coupling constants from 
the two approaches. 



VI. CONCLUSIONS 



We performed a resummation of the Gross-Llewellyn Smith (GLS) sum rule by fully accounting for the correct 
known form of the leading infrared (IR) renormalon singularity at b = 1 in the Borel plane. As one direct conse- 
quence of this singularity, the resummed "perturbative" part Ap (a s (Q)) of the GLS has a branch cut of the form 

(-a s (Q 2 ))~^exp(-7r/A)a s (Q))[l + 0(a s )], i.e., a twist-four term (-l)"" {1/Q 2 )a] 2 > ' Po [1 + 0(a s )} with the branch 
cut discontinuity factor exp(±i7w). Here, 72 =8/9 is the known one-loop coefficient of the anomalous dimension of 
the corresponding twist-four d = 2 operator appearing in the OPE. A crucial element of the analysis was the fixing of 
the "nonperturbative" part as the negative of the aforementioned branch cut term, thus making the resummed ( "per- 
turbative" plus "nonperturbative" ) GLS sum rule manifestly real. This procedure is free from the known ambiguity of 
separation of the "perturbative" and "nonperturbative" parts. In the Borel resummation of the "perturbative" part, 
we further employed a conformal transformation to minimize the numerical influence of other renormalon singularities. 
All this allowed us to perform the resummation of the massless part of the GLS sum rule, i.e., of the contributions 
of the three light quark flavors. The contributions of the heavy (c-)quark were added separately, as were the target 
nuclear correction contributions and the light-by-light contributions. These three types of contributions turned out 
to be small, in contrast to the "nonperturbative" contributions which turned out to be significant. The calculations 
were performed in the MS renormalization scheme, and the renormalization scale /z 2 (~ Q 2 ) was taken in the region 
of the smallest /i 2 -sensitivity of our results. 

We then confronted the resummed expressions with the Fermilab CCFR Collaboration data (l| for the GLS (at 
Q 2 = 2,3.16 GeV 2 ) which already include the target mass corrections. Our central value prediction for the QCD 
coupling parameter, corresponding to the central GLS values of the CCFR, is a s (Mz) ~ 0.117 [see Eqs. (|38| ) and 
([4l|)], in good agreement with the present world average. This is different from the central value predictions of 
previous analyses of the GLS sum rule by the CCFR Collaboration |l^] (a s (Mz) ~ 0.114) and by the authors of f46j ] 
(cts(Mz) ~ 0.115) which are below or at the lower edge of the world average. We have seen that our approach to the 
calculation of the nonperturbative contribution is consistent with the OPE approach, and the main reason for the 
difference between the two approaches is that at small Q 2 < 4 GeV 2 , the OPE approach tends to overestimate and 
the Borel resummed perturbative contribution cannot be approximated by an OPE amplitude. 

The GLS sum rule, at the low gauge boson transfer momenta Q 2 — 2-4 GeV 2 , is a very important quantity to 
measure, because it has apparently a strong nonperturbative component, stronger than in some other low-energy 
QCD observables such as the semihadronic tau decay rate. The more precise experimental values of the GLS sum 
rule would help determining the higher-twist contributions more accurately. 
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Appendix A. BRANCH CUT SINGULARITY AND THE NONPERTURBATIVE PART 



In this Appendix, we show explicitly formulas leading to Eqs. ([l0|)-( 12), using the identity (^|) and the leading IR 
renormalon singularity structure (0) . The latter structure around b = 1 can be rewritten more explicitly as 



Wo (b) = 



C 



Pl +V (\ - by 1 -" [1 + « a (l-6) + k 2 (1 - bf + ■ ■ ■] 



+ (Analytic part) 



(A.l) 
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The leading part is known (y), while the coefficients Kj (j > 1) of the subleading parts are not yet known. Inserting 
the above expansion into the Borel integration formula ((^) and performing the change of the integration variable 
b = 1 + float, we obtain 



ImAp (a(Q) ± ie) 



C 



r(-iz) 



where 



/>oo±ie 

/±(a) = Im/ dt e-\-t)- l - v [ 



-l//3oa(Q) a(Qr , /±(a(g)) 



1 + Ki(/3 a)(-t) + K 2 ((3oa) 2 (-t) 2 + ■ 



(A.2) 



= =F sin(7ri/) / 
Jo 



dt e~H~ 



[1 + «i(fl)0)(-t) + K 2 (/9 a) 2 (-t) 2 + 



= Tr(-y) sin(7rz/) [l + Ki(/3 a)f + « 2 (/3oa) M 1 ' -1 ) + ' 



(A.3) 



Now requiring the imaginary parts of Ap at positive a(Q) to match those in ( A.2 ) Ap for general complex a(Q) can 
be written as 

A P (o(Q)) = -Ce- 1 /^(Q)| ( „ a(g)) -, + Kl{ _ v)M _ a{Q)) -u+i 

+A P (a(Q)) , (A.4) 

where Ap(a) is the part with no singularities (no cuts) for a > 0. According to Ref. [^BJ, the nonperturbative 
part Anp(o(Q) ± ie) must cancel the imaginary part of Ap(a(Q) ± ie), and Anp wa s ch osen to have the simplest, 
presumably the most natural, form - i.e., just the negative of the branch cut term of ( A.4 ) 



(o(Q)) - +Ce- 1 /^ a W{(~a(Q))-» + Kl (- v )(3 (- a (Q))-»+ l 



(A.5) 



Further, since (— a^Fie) v+n = a v+n exp[±i(z/— n)7r], expressions ( ]A.2| ) and (A.5) immediately relate ReANp(a±ie) 
with the (calculable) ImAp (a ± ie) 



ReANp(a ± ie) = =p cot(i/7r)ImAp(a ± ie) , 



(A.6) 



giving the result (|l§). 



17 The absolute value of the singular term in (A.4) can be rewritten as ~ (A 2 /Q 2 )a(Q) +J2/l3o [l + O(a)], where A is a Q- 
independent scale, and 72 is the one-loop coefficient of the anomalous dimension of the twist-four (d=2) operator <C 0(Q) 3>~- 



A 2 a(Q) +72//3 ° appearing in the OPE 



15 



F. David, Nucl. Phys. B234, 237 (1984). 

F. David, Nucl. Phys. B263, 637 (1986). 

A. H. Mueller, Phys. Lett. B308, 355 (1993). 

G. Grunberg, Phys. Lett. B325, 441 (1994). 

D. E. Soper and L. R. Surguladze, Phys. Rev. D54, 4566 (1996), hep-ph/9511258. 
T. Lee, Phys. Rev. D56, 1091 (1997), hep-th/9611010. 

J. C. Le Guillou and J. Zinn- Justin, Large order behavior of perturbation theory, (Amsterdam, Netherlands: North- Holland, 
1990), and references therein. 

A. H. Mueller, The QCD perturbation series. Talk given at Workshop on QCD: 20 Years Later, Aachen, Germany, 9-13 
Jun 1992:162-171 (QCD161:W586:1992) . 

G. Cvetic and T. Lee, Phys. Rev. D64, 014030 (2001), hep-ph/0101297. 

I. Caprini and J. Fischer, Phys. Rev. D60, 054014 (1999), hep-ph/9811367. 

G. Cvetic, C. Dib, T. Lee, and I. Schmidt, Phys. Rev. D64, 093016 (2001), hep-ph/0106024. 

G. Cvetic, T. Lee, and I. Schmidt, Phys. Lett. B520, 222 (2001), hep-ph/0107069. 

J. H. Kim et ai, Phys. Rev. Lett. 81, 3595 (1998), hep-ex/9808015. 

J. H. Kim, "A measurement of a s from the Gross-Llewellyn Smith sum rule," PhD thesis, Columbia University, 1998 
(unpublished). 

T. Lee (2001), hep-ph/0104306, to appear in Phys. Rev. D. 

E. V. Shuryak and A. I. Vainshtein, Nucl. Phys. B 199, 451 (1982); ibid. 201, 141 (1982). 
V.M. Braun, and A.V. Kolesnichenko, Nucl. Phys.. B283, 723 (1987). 

B. Lautrup, Phys. Lett. B69, 109 (1977). 

V.A. Novikov, M.A. Shifman, A.I. Vainshtein, and V.I. Zakharov, Nucl. Phys. B249, 445 (1985). 

A comprehensive review on power corrections and renormalons can be found in: M. Beneke and V.M. Braun, hep- 
ph/0010208. 

A. L. Kataev, A. V. Kotikov, G. Parente and A. V. Sidorov, Phys. Lett. B 41 7, 374 (1998) |arXiv:hep-ph/9706534[ ; 

A. L. Kataev, G. Parente an d A. V. Sidorov, Nucl Phys. B 573, 405 (2000) [ |arXiv:hcp-ph/9905310| ]; A. L. Kataev, 
G. Parente and A. V. Sidorov, |arXiv:hep-ph/010622l| . 

U.K. Yang, and A. Bodek, Phys. Rev. Lett. 82, 2467 (1999); U.K. Yang, A. Bodek, and Q. Fan, |hep-ph/9806457 . 

B. V. Geshkenbein, B.L. Ioffe, and K.N. Zyablyuk, Phys. Rev. D64, 093009 (2001). 
A. Shifman, A.I. Vainshtein, and V.I. Zakharov, Nucl. Phys. B147, 448 (1979). 

T. Lee, Phys. Lett. B462, 1 (1999). 

G. Parisi, Nucl. Phys. B150, 163 (1979). 

C. N. Lovett-Turner and C. J. Maxwell, Nucl. Phys. B452, 188 (1995). 
M. Beneke, Phys. Rept. 317, 1 (1999), and references therein. 

O. Nachtmann, Nucl. Phys. B63, 237 (1973). 
P. M. Stevenson, Phys. Rev. D23, 2916 (1981). 

S. G. Gorishny and S. A. Larin, Phys. Lett. B 172, 109 (1986); E. B. Zijlstra and W. Van Neerven, Phys. Lett. B 297, 
377 (1992). 

5. A. Larin and J. A. M Vermaseren, P hys. Lett. B 259, 345 (1991); 
K.S. Jeong and T. Lee, |hep-ph/020415C . 

G. Grunberg, Phys. Lett. 95B, 70 (1980), HOB, 501(E) (1982); 114B, 271 (1982); Phys. Rev. D 29, 2315 (1984). 

A. L. Kataev, N V. Krasnikov and A. A. Pivovarov, Nucl. Phys. B 198, 508 (1982) [Erratum-ibid. B 490, 505 (1982)] 

[ bep-ph/9612326^ . 

A. Dhar and V. Gupta, Phys. Rev. D 29, 2822 (1984); V. Gupta, D. V. Shirkov, and O. V. Tarasov, Int. J. Mod. Phys. A 

6, 3381 (1991). 

A. L. Kataev and V. V. Starshenko, Mod. Phys. Lett. A 10, 235 (1995). 

G. Cvetic, Phys. Lett. B 486, 100 (2000); G. Cvetic and R. Kogerler, Phys. Rev. D63 (2001) 056013. 
S. J. Brodsky, G. P. Lep age and P. B. Mackenzi e, Phys. Rev. D 28, 228 (1983). A. L. Kataev and V. V. Starshenko, Phys. 
Rev. D 52, 402 (1995) |arXiv:hep-ph/9412305l; J. R. Ellis, I. Jack, D. R. Jones, M. Karliner and M. A. Sa muel, Phys 
Rev. D 57 2665 (1998) larXiv:hep-ph/9710302[ . F. A. Chishtie and V. Elias, Phys. Lett. B 521, 434 (2001) |arXiv:hep- 



ph/0107052 1. 

[40] M. Buza, Y. Matiounine, J. Smith and W. L. van Neerven, Nucl. P hys. B 485, 420 (1997) JarXiv:hep-ph/9608342 |. 
[41] J. Blumlein and W. L. van Neerven, Phys. Lett. B 450, 417 (1999) [ |arXiv:hep-ph/981135l| ]. 
[42] S. A. Kulagin, Nucl. Phys. A640, 435 (1998). 

[43] Particle Data Group, D. E. Groom et ai, Eur. Phys. J. C 15, 695 (2000). 



16 



T. van Ritbergen, J. A. M. Vermaseren and S. A. Larin, Phys. Lett. B 400, 379 (1997). 
K. G. Chetyrkin, B. A. Kniehl, and M. Steinhauser, Phys. Rev. Lett. 79, 2184 (19 97). 
J. Chyla and A. L. Kataev, Phys. Lett. B 297, 385 (1992) |arXiv:hep-ph/920921S ] . 
W. J. Marciano, Phys. Rev. D 29, 580 (1984). 

M. Dasgupta and B. R. Webber, Phys. Lett. B 382, 273 (1996) ]arXiv:hep-ph/960438cfl . 

Y. L. Dokshitzer, G. Marchesini and B. R. Webber, Nucl. Phys. B 469, 93 (1996) |arXiv:hep-ph/951233e[ ] 
S. Fajfer and R. J. Oakes, Phys. Lett. B 163, 385 (1985). 

I. Hinchliffe and A. V. Manohar, Ann. Rev. Nucl. Part. Sci. 5 0, 643 (2000) ]arXiv:hep-ph/000418r| . 
S. Bethke, J. Phys. G 26, R27 (2000) |arXiv:hep-ex/000402lt 
G. G 



Ross and R. G. Roberts, Phys. Lett. B 322, 425 (1994) |arXiv:hep-ph/9312237 




FIG. 1. Integration in the w— plane along the C_ ray w = x exp(— \<j>) (0<x<l, <j> = 0.67967) gives the same result as the integration 
parallel to the positive real axis (0 < to < 1) and arc w = exp(—i(j>') (0 < (j>' < </>). If integrating in the first quadrant, the paths are those 
obtained from the presented paths by reflection across the real axis, and the C-|_ ray is w = x exp(+i</>) ■ 
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FIG. 2. The (total) A(Q 2 ) = A(Q 2 ,P + NP) + A(Q 2 ,/./.) + A(Q 2 ,c - quark), as well as the n f = 3 perturbative part A(Q 2 , P) and 
nonpcrturbative part A(Q 2 ,NP), as functions of the renormalization scale /i 2 , as given by the applied resummation method. Given are 
the results at the N 3 L level, and for comparison, at the NNL level. The curves are for Q 2 = 2.00 GeV 2 and a s (Q 2 ,MS) = 0.3483. 
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FIG. 3. Same as in Fig. 0, this time for Q 2 = 3.16 GeV 2 and a a (Q 2 , MS) = 0.3053. 
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FIG. 4. The (total) A(Q 2 ) = A(Q 2 , P + NP) + A(Q 2 , 1. 1.) + A(Q 2 ,c - quark), and the separate parts, as functions of a s (Q 2 , MS), as 
given by the applied resummation method. The renormalization scale was fixed to be £ 2 = fi 2 /Q 2 = 3.3, and the VF-boson virtuality is 
Q 2 = 2.00 GeV 2 . The present experimental bounds and the central value, for the total A(Q 2 ) in this case, are denoted as three horizontal 
dashed lines. 
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FIG. 6. Borel resummed and OPE amplitudes versus Q 2 (GeV 2 ). ct s (2GcV 2 ) = 0.35 is assumed, (a): Borel resummed amplitudes 
of the perturbativc part Rc[A P (Q 2 )] at NLO, NNLO, and N 3 LO; (b): NLO (dot-dashed), NNLO (dashed), and N 3 LO TPS (solid) of 
Wo(a s (Q)); (c): Borel resummed A P (Q 2 )+A N p(Q 2 ) (solid line) against NLO (dashed) and NNLO (dot-dashed) OPE amplitudes. Dotted 
line denotes the Borel resummed with the wrong sign; (d): Borel resummed Re[Ap(Q 2 )] (solid) versus an NLO OPE amplitude (dashed). 
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FIG. 7. (a): N 3 LO Borcl transforms. The solid line represents the real part of the Borel transform with the rcnormalon at b = 1 
properly taken into account, and the dashed line represents the TPS Borel transform, (b): N 3 LO Borel resummed Re[Ap(Q 2 )] (solid line) 
versus N 3 LO TPS (dashed). 
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